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Abstract. Galileo, in the XVII century, observed that the small os- 
cillations of a pendulum seem to have constant period. In fact, the 
Taylor expansion of the period map of the pendulum is constant up 
to second order in the initial angular velocity around the stable equi- 
librium. It is well known that, for small oscillations of the pendulum 
and small intervals of time, the dynamics of the pendulum can be ap- 
proximated by the dynamics of the harmonic oscillator. We study the 
dynamics of a family of mechanical systems that includes the pendulum 
at small neighbourhoods of the equilibrium but after long intervals of 
time so that the second order term of the period map can no longer 
be neglected. We analyze such dynamical behaviour through a renor- 
malization scheme acting on the dynamics of this family of mechanical 
systems. The main theorem states that the asymptotic limit of this 
renormalization scheme is universal: it is the same for all the elements 
in the considered class of mechanical systems. As a consequence, we 
obtain an universal asymptotic focal decomposition for this family of 
mechanical systems. This paper is intended to be the first of a series of 
articles aiming at a semiclassical quantization of systems of the pendu- 
lum type as a natural application of the focal decomposition associated 
to the two-point boundary value problem. 



1. Introduction 

1.1. Focal decomposition. The concept of focal decomposition was intro- 
duced by Peixoto in [Ij (under the name of u-decomposition) , and further 
developed by Peixoto and Thom in [2]. The starting point was the 2-point 
boundary value problem for ordinary differential equations of the second 
order 

X = f{t, x,x) ,t, X, X £ M. 

x{ti) = Xi , x{t2) = X2 , (1) 

which was formulated for the equation of Euler at the beginnings of the 
calculus of variations in the first half of the 18th Century. This is the 
simplest and oldest of all boundary value problems. Accordingly, there is a 
vast literature about it, mostly in the context of applied mathematics where 
one frequently uses the methods of functional analysis. Here we adopt a 
different point of view: we look for the number i € {0, 1, 2, oo} = N of 
solutions of problem ([1]), and how this number varies with the endpoints 
{ti,xi) and {t2,X2). 
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Let = x M?{t2,X2) be the set of all pah's of points of the 

(i,a;)-plane and to each point {ti, xi,t2, X2) G associate the number of 
solutions i of the boundary value problem ([T|). When ti = t2, the index i 
is defined as being if xi 7^ X2, and 00 if j;i = X2- The diagonals 5, A are 
defined as 

6 = {{h,xi,t2,X2) eR^ :h = t2} 

A = {{ti,Xi,t2,X2) -.ti = t2, Xi = X2} . 

Let Sj C be the set of points to which the index i has been assigned. 
Clearly is the disjoint union of all the sets Sj, that is, these Sj define 
a partition of M^. This partition is called the focal decomposition of 
associated with the boundary value problem ([T|): 

M^ = i;oUSiu...uSoo . 

The decomposition above is the central object of study for a given /. 

From the above definition, the diagonal 5 is 3-dimensional and is naturally 
decomposed into the 2-dimensional connected component A belonging to 
and two 3-dimensional connected components belonging to Sq. This is 
called the natural stratification of 6. 

If one of the endpoints in ([l]) is kept fixed, say {ti,xi) = (0,0), then 
the sets Sj induce a decomposition of M^(t2,a^2) by the sets fij = Sj n 
({(0, 0)} X M^(t2, 3^2)) • The restricted problem with base point (0, 0) consists 
of finding the corresponding focal decomposition of by the sets af. 

M2 = O-Q U (Ti U ... U (Too ■ 

Throughout this paper, when we speak of the focal decomposition associated 
to problem ([T]), we will be referring to the focal decomposition of defined 
above. 

One has some clarification of what the sets fij, Sj are if one relates them 
to a certain surface associated to each point {t, x) called the star of this 
point. It is the lifted manifold introduced in [3j for other purposes. Let 

dt _ dx du „, , 

^ = 1 ,-r = u ,-r = fit, X, u) (2) 
dt dt ' dt •'^ ' ' ^ ' 

be the first system in R?{t,x,u) equivalent to ([T]). One calls the star asso- 
ciated to the base point (t, x) the surface S{t,x) C W^{t,x,u) obtained by 
the union of the trajectories of ^ passing through the points of the line 
{t,x,u), —00 < u < 00. One verifies immediately that if vr is the projection 
7r{t,x,u) = {t,x) then: a) with respect to base point {t2,X2) G o"i 

iff [vr I S{ti, xi)]~^ {t2, X2) consists of i points; also b) (ti, xi,t2, X2) G Sj iff 
S{ti,xi)riS{t2,X2) consists of i solutions of See [21/4] for more details on 
the relationship between the focal decomposition and stratification theory. 

Concerning the star, a final comment seems appropriate. To each curve 
7 solution of a two point boundary value problem, we associate two stars, 
one at each endpoint, which do intersect. When they do so transversally, 
we say that 7 is transversal or "structurally stable" because 7 can not be 
perturbed away by small variations of the equation. Besides, this situation 
is generic [3j, i.e. any equation can be approximated by one for which the 
solutions of the 2-point boundary problem are "structurally stable" . 
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A notorious example of a focal decomposition due to Peixoto and Thom 
[2], is provided by the focal decomposition of the pendulum equation x + 
sinx = with base point (0,0) (see Figured]). This focal decomposition 
contains non-empty sets cjj with all finite indices. Every set o"2fc-i, k = 
1,2,..., consists of a 2-dimensional open set plus the cusp-point (±/c7r,0); 
they all have two connected components. All four connected components 
of the even-indexed sets (T2k are open-arcs, asymptotic to one of the lines 
X = zb-TT and incident to the cusp-points (ibA;7r, 0); the lines x = ivr are part 
of (Ti, except for the points (0, itvr) which belong to do. 
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Figure 1: The pendulum's focal decomposition. 



There are many examples in the calculus of variations where rules are 
given for determining the number of solutions of the boundary value problem 
([T]) and for which the corresponding focal decomposition can be derived. The 
first such example seems to be a theorem of Goldschmidt (1831) concerning 
the number of catenaries passing through two points in the plane [Sj p. 
340]. See also the papers by MacNeish [B] and Peixoto [3], and especially 
the textbook by Collatz [7, pp. 170ff.], where several focal decompositions 
are given. However, the number of distinct indices obtained for the focal 
decompositions in these examples is very small and the concept itself is not 
recognized. 

After its recognition and early development in [U [2] , the concept of focal 
decomposition was further developed by Kupka and Peixoto in ^ in the con- 
text of geodesies and by Peixoto and Silva in [4] regarding the Euler equation 
associated with a given Lagrangian function L. In [9], Kupka, Peixoto and 
Pugh proved generic results pointing out the existence of residual sets in 
the space of complete Riemannian metrics on a finite dimensional complete 
Riemannian manifold M, for which bounds upon the number of geodesies 



4 C. A. A. DE CARVALHO, M. M. PEIXOTO, D. PINHEIRO, AND A. A. PINTO 

of prescribed length on M were obtained. Relationships with the arithmetic 
of positive definite quadratic forms are considered by Kupka and Peixoto in 
[8] and by Peixoto in [lOl [11] , and in [HI [10] attention is also drawn to the 
relationship between focal decomposition and the Brillouin zones of solid 
state physics, a view that is further developed later by Veerman et al. in 
[12] . In [11], it is pointed out that the focal decomposition associated to 
([1|) is relevant for the computation of the semiclassical quantization of this 
equation via the Feynman path integral method and in |13l [H] de Carvalho 
et al. exhibit further relations with quantum statistical mechanics. 

The concept of focal decomposition might also be relevant to the study 
of caustic formation by focusing wavefronts. This appears in optics (see the 
paper by Berry and Upstill [15] ) , tsunami formation (see the papers by Berry 
[TBI [T7] ) or general relativity (see the papers by Friedrich and Stewart [18] . 
Hasse et al. [E], Ellis et al. [20] and Ehlers and Newman plj). It is not 
so surprising that focal decomposition, i.e. the two-point boundary- value 
problem, is intimately connected to the very formulation of a variational 
problem, as Euler taught us. 

1.2. Renormalization. Besides focal decomposition, the other key topic 
in this paper is renormalization — the study of asymptotic self-similarity. 
The main idea behind renormalization is the introduction of an operator — 
the renormalization operator — on a space of systems whose action on each 
system is to remove its small scale behaviour and to rescale the remaining 
variables to preserve some normalization. If a system converges to some 
limiting behaviour under iteration of the renormalization operator then we 
say that such behaviour is universal. Since the renormalization operator 
relates different scales, such universal behaviour is self-similar. 

In this paper we introduce a new kind of renormalization scheme acting on 
a family of mechanical systems which includes the pendulum. This scheme 
has the distinguishing feature that time is not rescaled, hut rather translated, 
while the initial velocities and space are appropriately scaled. 

Renormalization is extremely relevant in several areas of physics. The 
first ideas concerning renormalization were introduced in the 1940's in the 
context of quantum electrodynamics by Bethe, Feynman, Schwinger and 
Dyson. These ideas were later developed by Stueckelberg and Petermann in 
[22] and by Gell-Mann and Low in [23] , with the introduction of the Renor- 
malization Group as a tool to improve approximate solutions to quantum 
field theory equations. Later developments in the subject were due to Callan 
[24| . Symanzik [25] and Weinberg [2B], among many others. Kadanoff in [27] 
and Wilson in [^ introduced renormalization techniques into statistical me- 
chanics in order to improve the understanding of critical phenomena. 

The relevance of renormalization in dynamical systems was first noticed 
by Feigenbaum in [23 [30] and, independently, by Coullet and Tresser in 
[31| . with their discovery of period doubling universality in unimodal maps 
of the interval. They introduced a renormalization operator — the period 
doubling operator — to show that period doubling sequences for this class 
of maps are asymptotically self-similar and that these sequences have an 
identical form for a large open set of such class of maps. This renormaliza- 
tion operator consists of time and space asymptotic normalizations chosen 
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to preserve the dynamical characteristics of the maps under iteration. Infin- 
itely renormalizable maps separate the regions with regular (zero entropy) 
and chaotic (positive entropy) dynamics and, moreover, their invariant set 
is universal, i.e. it is the same for all infinitely renormalizable unimodal 
maps. It was soon realized that the period-doubling operator was just a 
restriction of another operator acting on the space of unimodal maps ~ the 
renormalization operator - whose dynamical behaviour is much richer. 

Renormalization has been a very active research area in dynamical sys- 
tems in the past two decades: Sullivan in [32l [33], McMullen in [3l], de 
Melo and Pinto in [35], Lyubich in [36] and Faria, de Melo and Pinto in [37] 
studied the renormalization operator acting in families of unimodal maps. 
The renormalization operator naturally appears in several other families of 
maps, such as families of critical circle maps or families of annulus maps 
(see, for instance, MacKay [MIES], Ostlund et al. [10], Lanford [HI [12], 
de Melo [33], Martens [H], de Faria and de Melo [Ml [16] and Yampolsky 

[13 is]). 

The main subject of this paper is a renormalization scheme acting on 
the dynamics of a family of mechanical systems that includes the pendu- 
lum. Our motivation for the introduction of such a scheme comes from the 
restricted focal decomposition with base point (0, 0) of the pendulum equa- 
tion X + sin X = in Figure [TJ It turns out that the sequence formed by the 
even-indexed sets in the pendulum's focal decomposition is approximately 
self-similar. The renormalization scheme we introduce can then be justified 
in the following way: for a large integer n, we consider the even-indexed 
set (J2n and, contrary to previous renormalizations, we do not rescale time 
but just shift it so that its origin is at f = nvr. We then restrict the initial 
velocities to a small interval so that that the index corresponding to the 
shifted even-indexed set is equal to one; we complete the procedure by nor- 
malizing space in such way a that the shifted even-indexed set is asymptotic 
to the lines x = ±1. Under iteration of this renormalization scheme, we ob- 
tain asymptotic trajectories that define an asymptotic focal decomposition. 
Both the asymptotic trajectories and focal decomposition are universal and 
self-similar. 

To be more precise, our renormalization scheme acts on the dynamics of a 
family of mechanical systems (see the books |49 tl50l[5T| l52j for details on me- 
chanical systems) defined by a Lagrangian function L{x,x) = — V{x), 
where V{x) is a given non-isochronous potential, i.e. not all the periodic 
solutions of the corresponding Euler-Lagrange equation have the same pe- 
riod. For more information on isochronous potentials, see the paper [53] by 
Bolotin and MacKay and references therein. We obtain that the asymp- 
totic limit of this renormalization scheme is universal and as a consequence 
we obtain an universal asymptotic focal decomposition for this family of 
mechanical systems. 

This paper is the first step towards a broader research program, proposed 
by Peixoto and Pinto, connecting renormalization techniques, focal decom- 
position of differential equations and semiclassical physics. The next steps 
of this research program include a prove of the convergence of renormalized 
focal decompositions to the universal asymptotic focal decomposition and an 
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extension of this renormalization scheme to obtain a 4-dimensional universal 
asymptotic focal decomposition, i.e. with no restrictions on the base point 
of the boundary value problem. The ultimate goal of such program would be 
to deal with applications of this renormalization procedure to semiclassical 
physics. 

The present paper is structured in the following way: in section [21 we 
show why we believe that our program is relevant for semiclassical physics; 
in section [3l we introduce the main concepts we deal with throughout the 
paper and rigorously state our main results, while section S] is devoted to 
providing the reader with a global picture of the strategy of the proof of the 
main results. The remaining sections are more technical: section [5] deals with 
mechanical systems defined by some quartic potentials - we compute explicit 
solutions of the Euler-Lagrange equation associated with such mechanical 
systems in terms of Jacobian elliptic functions (see the textbooks [5 H l55 | l56] 
for further details on Jacobian elliptic functions) and use known estimates 
for such functions to obtain estimates for our explicit solutions; section [6] 
is devoted to the study of mechanical systems defined by perturbations of 
the quartic potentials - we use an implicit function argument inspired in a 
paper by Bishnani and MacKay [ST] to obtain estimates for the solutions of 
this second class of mechanical systems as a continuation of the estimates 
obtained previously for the explicit solutions associated with the quartic 
potentials; section[7]is devoted to the study of the period map for the periodic 
orbits of this family of mechanical systems; sections [H] and are devoted to 
the renormalization scheme introduced in this paper. We summarize in 
section [lOl 



Focal decomposition is in fact a first step towards semiclassical quan- 
tization. This was already recognized in the semiclassical calculation of 
partition functions for quantum mechanical systems, where the need to con- 
sider a varying number of classical solutions in different temperature regimes 
became evident [58] . 

We shall illustrate this for two physical quantities which can be expressed 
in terms of Feynman path integrals: in quantum mechanics, we shall con- 
sider the one-dimensional propagator between an initial point xi, at time ti, 
and a final point X2, at time t2', in quantum statistical mechanics, we shall 
study the one-dimensional thermal density matrix element in position repre- 
sentation. For further details on semiclassical quantization, see the textbook 
[59j by Feynman and Hibbs or the review paper [00] by DeWitt-Morette. 

In the quantum mechanical case, the propagator is just the time evolution 
operator exp[—iH{t2—ti)/h] computed between position eigenstates and 



2. Semiclassical physics 




G{xi,ti;x2,t2) = (X2I exp[-iii'(t2 -ti)/h]\xi) , 



where H is the Hamiltonian 
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The propagator may be written as a path integral 

r-z{t2)=X2 / ^r^- 

J z(ti)=xi \ 



h 



(3) 



where S is the classical action of the mechanical problem 



S[z] 



t2 



1 .2 

—mz 
2 



V{z) 



dt . 



The integral over z in ([3]) stands for a sum over all trajectories that connect 
xi at ti to X2 at t2, an object whose mathematical characterization has led 
to much investigation over the years. 

We shall only be interested in the leading semiclassical expression for the 
propagator, which can be formally derived from ([Sj). The resulting expres- 



sion IS 



Gsc{xi,ti;X2,t2) = {'^^^[Sn^]] exp ( 

n=l ^ 



IT 



The i in the upper limit is the same as in the focal decomposition, and 
indicates that the approximation is restricted to all classical trajectories, i.e. 
all solutions of the classical equation of motion which satisfy the boundary 
conditions z{ti) = xi and z{t2) = X2- The Sn in the exponential stands for 
the value of the action of the n^^ classical trajectory, whereas Sn'' denotes 
the second functional derivative of the action with respect to z{t), computed 
at the n^^ classical trajectory. The inverse square root of the determinant 
of that operator, the so-called van Vleck determinant, accounts for the first 
quantum corrections in a semiclassical expansion. 

In the case of quantum statistical mechanics, the thermal density matrix 
element is just the Boltzmann operator exp[— computed between posi- 
tion eigenstates \xi) and \x2), i.e. p{x2,xi) = (x2| exp[— The path 
integral for this quantity is given by 



p{x2,xi) = / [Dz{t)] exp 

Jz{0)=xi 

with the so-called euclidean action defined as 



S[z] 



S[z] 



f3h 



1 .2 

—mz 
2 



+ V{z) 



dr 



There is a crucial difference between this formula and that for the classical 
action of the mechanical problem: the sign in front of the potential. In fact, 
in the path integral formulation of quantum statistical mechanics one is 
led to investigate the mechanical problem defined by minus the potential. 
With this in mind, one may proceed along the same lines as in quantum 
mechanics to obtain a semiclassical approximation to the thermal density 
matrix element. It is given by 



,(x2,xi) = ^ |det[S'^2)]| exp 



n=l 



Sn 

IT 



Another important difference with respect to quantum mechanics is that 
the sum runs only over those solutions of the (euclidean) equation of motion 
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satisfying the boundary conditions which are local minima of the euclidean 
action S[z], whereas in quantum mechanics all solutions, i.e. any extremum, 
must be taken into account. That is why we use i' as an upper limit of the 
sum, defined as the number of extrema that are minima. Clearly, i' < i. In 
reality, we need a refinement of the focal decomposition to tell us which of 
the solutions are local minima. 

Either in quantum mechanics or in quantum statistical mechanics, the 
semiclassical approximation has to sum over all, or part of, the classical 
paths satisfying fixed point boundary conditions. Given the pairs {xi,ti) 
and {x2,t2), or {xi,0) and {x2,l3h), the number and type of classical tra- 
jectories are the very ingredients which lead to a focal decomposition. It 
should, therefore, be no surprise that the focal decomposition can be viewed 
as the starting point for a semiclassical calculation. 

As for the renormalization procedure, it was introduced to study the be- 
havior of classical trajectories for very short space and very long time sepa- 
rations of the fixed endpoints. It maps those trajectories into n-renormalized 
ones, whose time separations are shifted by n half-periods, and whose space 
separations arc scaled up to values of order one. As will be shown in the 
sequel, this procedure converges to an asymptotic universal family of trajec- 
tories that have a well-defined and simple functional form, and which define 
an asymptotic universal focal decomposition self-similar to the original one. 

The natural question to pose is whether the combination of focal decompo- 
sition and renormalization can be used to calculate semiclassical expansions 
for propagators in the short space, long time separation of the endpoints, 
or analogously, for thermal density matrices for short space separation and 
low temperatures (long euclidean time ph is equivalent to low temperatures 
T = l/{kBl3)) by using the simple asymptotic forms alluded to in the pre- 
vious paragraph. 

The conjecture to be investigated in a forthcoming article is that this can 
be done in a relatively simple way, thanks to the simple form of the asymp- 
totes. This will bypass a much more difficult (if not impossible) calculation 
involving Jacobi's elliptic functions. Should our expectation be realized, we 
would obtain semiclassical estimates for both propagators and thermal den- 
sity matrices in the short space/long time or short space/low temperature 
limits. Expressing those quantities in terms of energy eigenfunctions and 
energy eigenvalues, we have 

oo 

G{xi,ti;x2,t2) = ^ i/;^{x2)il^m{xi) ex.-p{-iEm{t2 -ti)/fi) , 

m=l 

in the case of quantum mechanics, or 

oo 

p{x2,xi) = ^ tp^{x2)'ilJm{xi) exp{-(3Em) (4) 

m=l 

in the case of quantum statistical mechanics. 

Clearly, in the latter case, if we take j3 large, only the lowest energy Eg 
(the ground state) will contribute. Furthermore, the points xi and X2 are 
to be taken in the limit of short space separation. If we choose one of them 
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to be the origin, then dH) becomes 

p{x2,xi) K. p(0,0) « |V'oPexp(-/3£;o) ■ 

This means that the combination of focal decomposition and renormaUzation 
may lead us to a direct estimate of the ground state energy for a quantum 
mechanical system from the asymptotic forms obtained in this article. 



3. Main Theorems and Definitions 

We start this section by fixing notation and introducing basic definitions 
which will be used throughout the paper. We also state the main results to 
be proved in the following sections. 

3.1. Setting. We consider mechanical systems defined by a Lagrangian 
function £ : ^ R of the form 

where the potential function V : M ^ M is a non-isochronous potential. 
Furthermore, we assume that the potential V is a C map (k > 5) with a 
Taylor expansion at a point G M given by 



V"(q*) 



±0[\q 



2 ' ' 4! 

where V"{q*) > and V^'^Xq*) / 0. The Euler-La grange equation associ- 
ated with dS]) is 



d'^q dV, , 

dr^ aq 

and the corresponding Hamilton equations are given by 

dq p 

dr m 

dp dV, , 

d7 = -d^(^)- 



(6) 



(7) 



Therefore, the point q* is an elliptic equilibrium of ([6]) (or equivalently, 
(g*,0) is an elliptic equilibrium of ([7])) and thus, there is a 1-parameter 
family of periodic orbits covering a neighbourhood of the equilibrium point. 



3.2. Asymptotic universal behaviour for the trajectories. Since q* 
is an elliptic equilibrium of ([6]) there is e > such that for all initial velocity 
V € [—€, e] the solutions g(z^; r) of the Euler-Lagrange equation ([6]) with 
initial conditions q{i';0) = q* and dg/dr(i^;0) = v are periodic. Thus, the 
trajectories q : [— e, e] x M ^ M of ([6]) are well-defined by q{u; r) for all r G M 
and u £ [— e, e]. Furthermore, there exist a > small enough and > 1 
large enough such that, for every n > N, the n-renormalized trajectories 
Xn '■ [—1, 1] X [0, an] — > R are well-defined by 



xn{v;t) = {-ly 



n,t H- 



q rn,t IJLU] v; 



nvr 



It 
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where Tn^t is the {n,t)-scaling parameter 

, 1/2 

(8) 



n,t 



8t ^ 



37rn 

£ = ±1 depending on the sign of V^^\q*) and oj and fi are given by 

Note that and fi are the natural time and length scales for the dynamical 
system defined by ([6]). Furthermore, the variables v and t are dimensionless, 
as well as the (n, t)-scaling parameter F^,*. Therefore, the n-renormalized 
trajectories Xn{v;t) are dimensionless. 

Definition The asymptotic trajectories Xi : [—1,1] x M.^ M are defined 

by 

Xe{v;t) = V sin {it (w^ - l)) , 
where i = ±1 depending on the sign of V^^\q*). 

The rest of the paper is mainly devoted to prove the following result and 
some of its consequences. 

Theorem 3.1. There exists f3 > small enough, such that, for every < 
e < 1/3, we have that 

-Xf(?;;t)||cro([_i,i]x[o,/3ni/3-^],iR) < O (n"^'/^) , 
where I is the sign ofV^^\q*). 

3.3. Asymptotic universal focal decomposition. The asymptotic tra- 
jectories Xi(v;t) induce an asymptotic focal decomposition of the cylinder 
C = M.Q X [—1, 1] by the sets cjj whose elements are pairs (t, x) G C such that 
Xi^Vjt) = X has exactly i solutions v{t,x) € [—1,1], each distinct solution 
corresponding to an asymptotic trajectory connecting the points (0, 0) G C 
and {t, x) £ C. Therefore, for each i £ {0, 1, oo}, the set o", C C contains 
all points in (t, x) G C such that there exist exactly i asymptotic trajectories 
connecting (0,0) G C and {t,x) G C. The following result is a consequence 
of theorem 13.11 

Theorem 3.2. There exists an asymptotic universal focal decomposition 
for the Euler-Lagrange equation ([6]) induced by the asymptotic trajectories 
Xi{v-t). 

The asymptotic universal focal decomposition is shown in Figure [21 As in 
the case of the focal decomposition with base point (0, 0) of the pendulum 
equation x + sin x = (see Figure [T]) , the asymptotic universal focal decom- 
position also exhibits non-empty sets cjj with all finite indices. However, 
contrary to that focal decomposition, which gives a stratification for the 
whole M^, our asymptotic universal focal decomposition gives only a strat- 
ification of the half cylinder C = x [—1, 1]. There are two main reasons 
for this to happen which we pass to explain. Firstly, our renormalization 
scheme acts only on periodic orbits, neglecting the high-energy non-periodic 
orbits, which restrains Xi^{v;t) to the interval [—1,1]. Secondly, we have 
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defined the renormalization operator only for positive times. Noticing that 
the mechanical systems we renormalize have time-reversal symmetry, one 
can extend the asymptotic universal focal decomposition to the cylinder 
R X [—1,1] by a symmetry on the x axis. 




X = -1 

Figure 2: The asymptotic universal focal decomposition. 



Some comments about Figure [2] are appropriate here. 
For every k € N, the set a2k is a 2-dimensional open set with two con- 
nected components. The the odd-indexed sets cr2k-i are the union of two 
open arcs, asymptotic to one of the lines x = ±1 and incident to the cusp- 
point ({k — l)7r,0), and a line segment joining the cusp points {{k — l)vr, 0) 
and (/c7r,0); the lines x = ±1 are part of do- 
Thus, the even-indexed sets (T2fc are 2-dimensional manifolds, while the 
odd-indexed sets cr2k-i are not manifolds because they contain the cusp- 
points {kn, 0). 

On the other hand, if we decompose the odd-indexed sets cr2k-i into a 
cusp-point plus three 1-dimensional manifolds (two open arcs and one line 
segment), then we get a decomposition of the whole plane into a collection 
of disjoint connected manifolds. The above decomposition of C is an ex- 
ample of what is called a stratification of C, the strata being the disjoint 
connected manifolds into which C was decomposed. Hence cJi consists of 
two 1-dimensional strata, a2 consists of two 2-dimensional strata, con- 
sists of three 1-dimensional strata and one 0-dimensional strata and so on. 
To complete the picture, ctq consists of two 2-dimensional strata plus the 
X-axis minus the origin which belongs to aoo- 
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4. The renormalization procedure 

In this section we introduce an affine change of coordinates of space and 
time that enables us to map trajectories of the Lagrangian system ([5]) to 
trajectories of a dimensionless Lagrangian system, simphfying the study of 
the asymptotic properties of such trajectories. 

4.1. An afRne change of coordinates and the perturbed quartic po- 
tentials. In order to simphfy the proofs throughout the paper, and without 
loss of generality, we apply a change of coordinates of space and time to the 
Lagrangian ^ and corresponding Euler-Lagrange equation The new 
dimensionless coordinates (x, t) are defined by 

X = (q-q*) , t = uoT , (10) 

where oj and /i are as given in Q. We obtain a normalized Lagrangian 
function 

L(x,x) = - V{x) , (11) 
with potential function V{x) of the form 

V{x) = V{S)) + \x^ + y^ + f{x) , (12) 

where 

a) f{x) G O {\x\^) is a C"^ map {k > 5); 

b) ^ = 1 if V^^\q*) > and ^ = -1 if V^^\q*) < 0. 

Since adding a constant to the potential does not change the form of the 
associated Euler-Lagrange equation, for simplicity of notation and without 
loss of generality we rescale the total energy associated with the Lagrangian 
system (lllh so that that V{0) = 0. In the next definition, we introduce the 
quartic potentials and the perturbed quartic potentials 

Definition The quartic potentials : M ^ M are defined by 

V,ix) = ^x^ + ^x\ (13) 
The perturbed quartic potentials : M ^ M are defined by 

V{x) = ^x' + ^-x' + f{x) , (14) 

where i £ {-1, 1} and \ f^^{x)\ G O {\x\^-') for i G {0, . . . , 5}. 

In the new coordinates {x, t), the Euler-Lagrange equation associated with 
the normalized Lagrangian (jlip reduces to 

i = -{x + lx^ + f'{x)) , (15) 

where the dot denotes differentiation with respect to the variable t. The 
next Lemma describes the relation between the solutions associated with 
the Lagrangian ([5]) and the normalized Lagrangian (jlip . The result follows 
easily from the change of coordinates (jlOp . 



AN ASYMPTOTIC UNIVERSAL FOCAL DECOMPOSITION 



13 




Vi{x) 



V-i{x) 



6 



X 







X 



1 



1 



-1 



1 



(a) 



(b) 



Figure 3: The quartic potentials Vf(x). 



Lemma 4.1. The solution q{i']T) of Euler- Lagrange equation ([6]) with ini- 
tial conditions q{iy; 0) = q* and ^{i^', 0) = is related to the solution x{v; t) 
of the Euler-Lagrange equation (jlSp with initial conditions x{v; 0) = and 
x{v; 0) = V by 



where uj and fi are as given in Q. 

4.2. The renormalization procedure. Since x = is a local minimum 
of the potential V{x) given by (fT2|) there is VMax > such that for all 
< VMax the solutions x{v;t) of the Euler-Lagrange equation (fT5|) with 
initial conditions x{v;0) = and x{v;0) = v are periodic. Furthermore, 
there exist a > small enough and > 1 large enough such that, for 
every n > N, the n-renormalized trajectories x„ : [—1, 1] x [0, an] R are 
well-defined by 



where r„^f is the (n, t)-scaling parameter defined in ([8]). Using Lemma l4.lt 
we obtain that Theorem 13.11 follows as a consequence of the following result. 

Proposition 4.2. There exist K > and N > 1 such that for every n > N 
and every t £ [0,Kn], the n-renormalized trajectories Xn{v;t) converge to 
the asymptotic trajectories Xi{v;t). Furthermore, the following bound is 
satisfied 



for all V e [-1, 1]. 

The rest of this paper is devoted to prepare the proof of Proposition 14. 2^ 
which will be completed in section [9l 



In this section we will make use of Jacobian elliptic functions (see [531 
[55| [56] ) to compute explicit solutions, with suitable initial conditions, for 
trajectories of a given class of Lagrangian systems. Furthermore, we will 
compute the first terms of the Taylor expansion of such solutions in terms 




xn{v;t) = {-ir r-j 



X {Tn,t v;mT - £t) 




5. QuARTIC POTENTIALS 
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of their initial velocities. In the next section, we will use these Taylor ex- 
pansions to study the trajectories associated to the more general potentials 
described in section HI 

We will consider Lagrangian functions of the form 

Le{x,x) = - Ve{x) , 

where Ve{x) is a quartic potential introduced in (fT3]l . and compute explicit 
solutions of the corresponding Euler-Lagrange equation 

x = —X — £x^ 

x(0) = (16) 
x(0) = V . 

We denote such solutions by X(:{v] t), where i S { — 1, 1}. 

Recall that x = is a local minimum of Vi{x). Thus, there exists VMax > 
such that for all initial velocities v satisfying \v\ < VMax, the solutions of 
([16]) are periodic. 

Let £ S {—1,1} and let u > be small enough so that the solutions of 
are periodic. We define the {£,v)- coefficients ai and bi by 

1/2\V2 



■= aiiv) = + + 2iv^) 

bi:=bfXv) = [-£ + l{l + 2lv'^f''^Y'^ ■ (1'^) 
Lemma 5.1. Let i = 1. The solution xi{v;t) of (|16|) is given by 

xi{v;t) = Aisd{Xit;mi) , 

where: 

(i) sd(t; m) is a Jacobian elliptic function. 

(ii) the amplitude A\, the frequency \\, and the parameter mi are given 
by 

Ai := Ai{v) = sign{v)aibi/ (ai^ + bi'^f^^ 
Xi:=\i{v) = ((ai2 + 6i2)/2)i/2 ^^g^ 
mi := mi{v) = (ai^ + . 

(iii) ai and bi are as given in (fTl 



Proof. Noting that the solution of (I16p is also the solution of the first order 
differential equation 

h' + Vi{x) = lv' 

with initial condition xi{0) = 0, integrating the above equation we obtain 
that the solution xi{v;t) is implicitly determined by 

-1/2 



Jo U^'^'^^^^V = sign(?;)t/V2 . 

from [S3], we get 
(sd~^ (xi/Ai; mi)) = sign(u)t/ 



Using formula 17.4.51 from [S3], we get 
1 
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where Ai, Ai and mi are defined in (jlSp in terms of oi and bi given in 
(fT7|) . The lemma follows by solving the above equality with respect to xi 
to obtain t). □ 

The proof of the next lemma is analogous to the proof of lemma ET] (using 
formula 17.4.45 from [54j ) . We skip its proof. 

Lemma 5.2. Let I = —1. The solution x^i{v;t) of (I16p with \v\ < VMax is 
given by 

; t) = sn (A_ii; m_i) , 

where: 

(i) sn(t; m) is the Jacobian elliptic sinus. 

(ii) the amplitude A^i, the frequency A_i and the parameter m_i are 
given by 



A. 


.1 := A^i{v) 




sign(w)&_i 


A- 


-1 ■■=X-i{v) 




a_i/\/2 


m_ 


1 := m_i(?j) 




(b-i/a^if 


are 


as given in 







In the two next lemmas we provide estimates for the amplitude Ag, fre- 
quency \i and parameter mi. The estimates are obtained by finding the 
values of the successive derivatives of Ai, Xi and me when v = 0. 

Lemma 5.3. Let i = 1. There exists VMax > such that for all v G 
{-VMax,Vmax), wc havc that 

Ai ev-\v'^±0[\v\^) 



mi G ± O [v^] . 

Lemma 5.4. Let i = —1. There exists VMax > such that for all v G 
{-VMax, VMax), we have that 

A^i ev + ^v^ ±0{\v\^) 

A_i G 1 - ± O {v^) 

m_i G ^-u^ ± O [v^] . 

Lemma 5.5. Let ji = ^ and j-i = — | and define ee{v) = 1 + jev"^ , so that 
we have 

A, = ee{v) ± 0{v^) . 
There is vq > such that for all \v\ < vq, the following estimates hold 

sm{Xi{v)t) G sm{ee{v)t) ± 0{tv^) 
sm.{ei>{v)t) G sin(t) ±0{tv'^) . 



16 C. A. A. DE CARVALHO, M. M. PEIXOTO, D. PINHEIRO, AND A. A. PINTO 

Proof. We use the estimates for in Lemmas 15.31 and 15.41 the definition of 
ei and the trigonometry formulas for the sine and cosine of the sum of two 
angles. □ 

The next Lemma gives us two classical estimates on Jacobian elliptic 
functions (formulas 16.13.1 and 16.13.3 of |54j ) . 

Lemma 5.6. There is mo > small enough such that for all \m\ < tuq, the 
following estimates hold 

sn(i; m) G sin(t) — -m {t — sin(t) cos(t)) cos(t) ± O (m^) 

dn(t; m) G 1 — — msin^(t) lb O (m^) . 

In the next Lemma we rewrite the classical estimates of Lemma 15.61 in a 
form that will be useful for the proof of Lemma 15.81 

Lemma 5.7. Let ge{t) be given by 

gi{t) = t cos(t) + i (sin(t) cos^(t) + 2 sin(t)) . 

There is tuq > such that for all \m\ < mo, the following estimates hold 

sd{t; m) G sin(t) - ^m (c/i (t) - 4 sin {t))±0 (m^) (19) 

sn(f;m) G sin(t) - ^m(5r_i (t) + 2sin(t)) ± O (m^) . (20) 

Proof. For the first inequality, using the definition of the Jacobian elliptic 
function sd(t;m) (see [M]), we obtain 

sn(t; m) 

sd{t; m = ( ■ 21 

^ ^ dn(t;m) ^ ' 

Putting together (j2ip and Lemma \5.6\ we get that there is mo > such that 
for all |m| < mo, we get 

sn(t; m) 



sd(t; m) 



dn(t; m) 

sin(t) — m{t — sin(t) cos(t)) cos(t)/4 ± O (m^) 



(22) 



l-msin^(t)/2±0 (m^ 
By Taylor's series, we get 

1 1 

1 - -m sin^it) ± O (m^) J C 1 + -m sin2(t) ± O (m^) . (23) 

Putting together inclusions ()22p and (1230 and the definition of gi{t), we get 

sd(t; m) G sin(t) - ^m {gi (t) - 4 sin (t)) ±0 (m^) , 

which proves (jl9p . 

The second equality follows from Lemma [5.61 and the definition of g-i (t). 
We obtain that there is mg > such that for all |m| < mo, we get 

sn (t; m) G sin (t) - ^m {g^i (t) + 2 sin (t)) ± O (m^) , 
as required. □ 
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Lemma 5.8. There is vq > such that for every \v\ < vq, the solution 
x^{v;t) of (fT6]) satisfies the following estimate 

XI {v; t)ev sm{ee{v)t) - lv^ge{ee{v)t) ± 0{t\vf) ± O {\vf) , 



where ee is defined in Lemma 15.51 and gi is defined in Lemma 5. ' 
Ln particular, for every i G { — 1, 1}, we get 

Xi{v;t) evsm{t)±0{t\v\^)±0{\v\^)) . 

Proof. We separate the proof into two cases: i = 1 and i = —1. 

Case i = 1. For simplicity of notation, we will use A = Ai, X = Xi, 
m = mi , g = gi and e = ei . 

From Lemma l5.H we have that 

xi{v;t) = Asd{Xt;m) . (24) 
Combining (jl9p in Lemma |5 . 71 with ()24p . we obtain 

xi{v;t) e A (^sm{Xt) -^m,{g{Xt)- 4 sin {Xt))±0{m^)^ (25) 

Combining the estimate for m in Lemma 15.31 with (|25p . we obtain 

xi{v;t) e A(^sm{Xt)-^(^^v^±0{v^)^{g{Xt)-Asm{Xt))±0{v^ 

C A(^sm{Xt)-^v^{g{Xt)-ism{Xt))±0{v^)^ . (26) 
Putting together the estimate for A in Lemma 15.31 and (j26p , we get 

xi{v;t) e (^v-^v^±0{\vf)^ (^sm{Xt)-^v'^{g{Xt)-Asm{Xt))±0{i 

C vsm{Xt) -lv^g{Xt)±0 {\vf) . 
8 

Applying Lemma 15.51 to ([27|) , we get 

xi{v;t) e vsin{et) -lv^g{et) ±0 {t\vf) ±0 {\v\^) 

8 

C vsm{et)-lv^g{t)±0{t\v\^)±0{\v\^) , (28) 
8 

as required. 

Case i = —1. For simplicity of notation, we will use A = A = A_i, 
m = m_i, g = g^i and e = e_i. 
From Lemma 15.21 we have that 



x_i{v;t) = Asn{Xt;m) . (29) 
Combining (I20p in Lemma 15.71 with (I29p . we get 



-iiv;t) = A (^sin (Xt) -^m{g{Xt) + 2 sin {Xt))±0 {m"^)^ . 



(30) 
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Combining the estimate for m in Lemma 15.41 and (j30p . we obtain 

x-i{v;t) e A(^sm{Xt) -J (^^v^ ±0 {v^)^ {g{Xt) + 2 sin {Xt))±0{v^)^ 

C A(^sm{Xt)-^v^{g{Xt) + 2sm{Xt))±0{v^)^ . (31) 
Putting together the estimate for A in Lemma 15.41 and ()3ip , we get 
x_i {v; t) G (^v + ^v^ ±0{\vf)^ (^sm{Xt) - ^v^ {g (Xt) + 2 sin (At)) ± 0{i 

C vsm{Xt)-lv^g{Xt)±0{\vf) . 
8 

Using Lemma 15.51 in ()32p . we obtain that 

x-i{v;t) G vsm{et)--v^g{et)±0{t\vf)±0{\v\^) 

8 

C 't;sin(et) - ^7;35(t) ±0(t|7;|5) ±0 (|w|5) , (33) 
8 

as required. 

Finally, applying again Lemma 15.51 to inequalities ([28l) and p3|) , for £ G 
{-1,1}, we get 

X£ (v; t) G t; sin(t) ± ©(t]-;;!^) ± O (l^l^) . 

□ 

6. Perturbed quartic potentials V 

In this section, we will consider perturbed quartic potentials V{x) that 
are perturbations of the quartic potentials V£{x) near the elliptic fixed point 
(0, 0) of the first order differential system associated to the Euler-Lagrange 
equation (ITS]) . We will prove estimates for the solutions x{v; t) of the Euler- 
Lagrange equation (fT5|) associated to V{x) using the solutions Xi{v;t) asso- 
ciated to the quartic potentials Vi{x). 

Along this paper, we will denote by x{v; t) the solutions of the Euler- 
Lagrange equation (fTS]) with initial conditions x{v; 0) = and x{v; 0) = v. 

Let {[a, b],R) be the set of all maps e {[a, b],R) such that V'(O) = 
0. Let U : ([0, T], M) ^ C° ([0, T], M) x M be the linear operator given by 

C/^ = (^^/, + ^,^(o) 

Lemma 6.1. The linear operator U : C^([0,r],M) C° {[0,T],R) is in- 
vertible and its inverse : C° ([0, T], M) x M ^ ([0,r],M) is a linear 
operator with bounded norm. 

Proof. The existence and uniqueness of (/(t),f) is guaranteed by the 
theorem of existence and uniqueness of solutions of differential equations. 
For every / G C° ([0,T],M) and every v G M, let us show that 

ipvit) = ■usin(t) — /(s) sin(s)ds^ cos(t) + /(s) cos(s)ds^ sin(t) 
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is equal to U~^{f, v). The first derivative of ipvit) is given by 
ijy{t) = f cos(t) + /(s) sin(s)ds^ sin(t) + /(s) cos(s)ds^ cos(t) . 
Tlie second derivative of ip^ (t) is given by 

^^(t) = -vsm{t)+f{t)+(^J /(s)sin(s)ds^ cos{t)-(^J /(s) cos(s)ds^ sin(t) 

Hence, ipvi't) + ipvit) = fit), ipv{t) = and ipyiO) = v which imphes that 
U{i>y{t)) = if{t),v). Therefore, the inverse U"^ : {[0,T],R) xR ^ 
Cq ([0, r],M) of [/ is a hnear operator, and 



\\U if{t)^v)\\cmo,T],R) - (1 + 2^) 11-^ IICO([0,Th 

Therefore, the hnear operator has bounded norm. 



+ V 



□ 



Theorem 6.2. For every T > 0, there exist 6i > 0, (52 > 0, > 
and K > with the following properties: For every \v\ < 6i and every 
e G Cq {[—82,62] ,R), with ||e||(72([-52 ^2] k) ^ ordinary differential 

equation 

X + X + ix^ + e{x) = 

has a unique solution Xe{v; t) G ([0, T] , M) with x^{v; 0) = and Xe(v; 0) = 
V, where £ = ±1. Furthermore, we have that 

Sx 

-?7. , . , , e <K\\e{x,{v;t))\\co^[o^T],R) ■ (^4) 

Proof Let 9 : {[0,T],R) x {[-6,6], R) ([0,T],M) x M, be given 

by 

e {x{t),e{x)) = {x + x + lx^ + e(2;), 2;(0)) . 
The nonhnear operator is C^, with partial derivative ^ (^x{t),t(x)) 

■■Cl{[Q,T], 

C°([0,T],IR) X R given by 

V = fV- + + 3^2:V + ^e(2;)V', V'(O)') , 

ox {x{t),e{x)) V / 

and with partial derivative ^(^^(J^^ ^(j,)) : Cq {[—6, 5] ,R 
given by 

a = (a, 0) . 



CO([0,r],R) X 



c°([o,r],M) X 



6e {x{t),e{x)) 

By Lemma 16.1^ the linear operator U : Cg ([0,T] 
given by 

Ui;= (V; + V,V'(0)) 
is invertible and its inverse has bounded norm. Furthermore, we have that 



6X {x(t),e{x)) 



CO([0,T],]R)xI 



< II (3£xV + De{x)i;, O) ||f.o([o,T],R)xR 

< (3C(522 + 53)||V'llc2([0,T],R) • 
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Hence, there is ^2 > and 6^ > small enough so that 



U 



5e 

dX {x{t),€{x)) 



< \\u- 



Therefore,the linear operator ^ : Cg ([0, T], 

vertible and its inverse has bounded norm 

^1 

r-l||-l 



6& 

6x (x{t),e{x)) 



-1 



u 



CO([0,r],M) X 



5Q 

6x {x(t),e{x)) 



IS m- 



such that 



By Lemmas 15.11 and 15.21 there exists X£{v;t) £ Cq{[0,T], 
@ {xi{v;t),0) = 0. Hence, by the implicit function theorem and by the 

5x ' 



invertibility of the operator there exist 6i > and ^3 > small enough 



with the following property: for all 1^1 < 61 and for all ||e||(^2([_^2 fe] 
there exists a unique solution Xe{v;t) £ Cq{[0,T],W) such that 

e{x,{v;t),e{x)) = . 
Differentiating (f35|l with respect to e, we get 



< <5,^ 



(35) 



se Sx_ 

5x {x^{v;t),e(x)) Sc (e{x)) 

Hence, the operator ^ : ([-5,(5] ,1 



+ 



5e 

6e {xciv;t),e{x)) 







is given by 



5x 

6e {e(x)) 



6e 

6x {x,{v;t),e{x)) 



6e 

de {x^{v;t),e{x)) 



Therefore, there is K > such that 



5x 

6e {e{x)) 



a 



(72([0,T],: 



60 

5x (x4v;t),e{x)) 



60 

6e {x^{v;t),e{x)) 



a 



c2([o,T],: 



< 



< 



60 

6x (x4v;t),e{x)) 

||a(x,(t;;t))||^o([o,T] 



60 

6e {xc{v;t),e{x)) 



a 



(S'O([0,T],IR)xI 



n 



Lemma 6.3. Let V{x) be a perturbed quartic potential. For every T > 0, 
there exists vq > such that, for all < t < T and for all \v\ < vq, the 
solution x{v; t) of the Euler-Lagrange equation associated to the potential 
V{x) satisfies the following estimate 

\\x{v-t) - a;K^;*)ll(72([o,T],iR) < O {v^) , 

where xe{v;t) are the solutions of the Euler-Lagrange equation associated to 
the potential Ve{x) and i is equal to the sign of V^'^\0). 

Proof. Let us consider the equation 

x + x + ex^ + e{x) = , (36) 

Let f'{x) be as given by (jl4p . Using Theorem 16.21 there exists vq > small 
enough such that, for every < k < 1 and every 1^1 < vq, the equation ([36]) . 
with e(x) = kf'{x), has a unique solution Xk{v;t) with initial conditions 
2^fc(^;0) = and ifc(f;0) = v. Furthermore, the solution Xk{v;t) is periodic 
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and there exists Ci > such that \xk{v;t)\ < Ci\v\ for all < t < T (see 
Lemma m]) . By p^ . there exists C2 > such that 



Hence, using ()34p . there exists K > such that 



\x{v;t) - xe{v;t)\\(^,2,,o,T],R) < sup 
"■^ ' " ' 0<k<l 



C2([0,T], 



5e fc/'(a;)'' 

< ^ l|/ (^fe(^'*))|lcO([0,T],R) 

< ^ ||e (a:)||co([_c'^|^|_c'j|^|]^]8) 

< KC2V^ . 

□ 

Proposition 6.4. Let V{x) be a perturbed quartic potential. For every 
T > 0, there exists vq > such that, for all < t <T and for all \v\ < vq, 
the solution x{v; t) of the Euler-Lagrange equation associated to the potential 
V{x) satisfies the following estimate 

X {v; t) ev sm{eg{v)t) - ^v^ ge{ee{v)t) ± O {v"^) , 

o 

where Xi{v;t) are the solutions of the Euler-Lagrange equation associated to 
the potential Vi{x) and i is equal to the sign of V^^\0). In particular, we 
get 

x{v;t) e vsm{t) ±0 {\vf\t\ +v'^) . 
Proof. The proof follows from combining Lemmas 15.81 and 16.31 □ 

7. Period map T 

In this section, we are going to compute the first three terms of the Tay- 
lor's expansion of the period map T{v) which associates to each initial ve- 
locity / of a solution x{v;t), with x(^;;0) = 0, its smallest period T{v) 
in a small neighbourhood of 0. 

Since x = is a local minimum of the perturbed quartic potential, the 
solutions x{v;t) of the Euler-Lagrange equation ([TS]) . with initial conditions 
x{v; 0) = and x{v; 0) = v, are periodic for all small values of v. Hence, 
there exists vq > 0, small enough, such that, for all \v\ < vq, the extreme 
points x"^{v) < < x^^ {v) of a solution x{v;t) are well-defined. Since the 
energy of the system is conserved along its orbits and the initial energy is 
equal to ^u^, we obtain that x"^{v) and x^^ (v) are implicitly determined by 

V{x^^{v)) = y(x"(i')) = ^v^ . 

Similarly, we will denote by xj^{v) < < xf^{v) the extreme points of the 
solution Xi{v;t) of ([TO]) . 

Definition The period map T : {—vq,vq) — > M is defined by 

_ fX^Hv) /. X -1/2 

T{v) = V2 -v"^ -V{x) dx . 

Jx"^(v) \2 / 
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We note that the period of the solution x{0, t) is equal to 0, but we keep 
T(0) = 2tt for the period map T to be continuous and smooth. 

Lemma 7.1. Let V{x) be a a perturbed quartic potential and let Vi^x) be 
a quartic potential. There exists vq > small enough such that, for all 
\v\ < vq, we have that x^{v) = —xf^{v), x'^{v) e 0(1^1), and 

x^^^ (v) e xf (v) ± O {\vf) andx''\v)exT{v)±0{\vf) . 

Proof. Since V{x) G ix^iO (x^) and V{x^iv)) = y(x™(?;)) = ^v^, we get 
E ±0{\v\) and E ±0{\v\). Recall that Ve{x) = ^x^ + ^x^. The 

solutions of the equation {x({v)) = ^v"^ are given by 

1 /9 

xf{v) = -xT{v) = (^-£ + i{l + 2£v^f^^^ 
By the definition of the perturbed quartic potentials in (fn|) . we have that 
V{x) = + + 

g ^X^^^X^±0{\xf) 

C Ve{x)±0{\xf) . 

Therefore, we get 

Vi{x^{v)) E V{x^{v))±0{\x^{v)\^) 

C lv'±0{\xf) . 



Hence, we have that 

C xf{v)±0{\vf) . 



x^^iv) E { -i + e{l + 4i{-v^ ±0{\xf) 



1/2^ 



1/2 



Similarly, we get 

x'^iv) E - (^-e + e(^i + A£(^^v'^ ±oi\xf] 

5\ 



1/2' 



1/2 



C xf{v)±0{\vf) . 



□ 



Lemma 7.2. Let V£{x) be a quartic potential. The period map T : {—vq, vq) 
M satisfies the following estimate 

Tiiv) = 27r - ^iv"^ ± O (v^) , 

where i is equal to the sign of V^'^\0). 
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Proof. By Lemma |7. 11 we have that 

.xf{v) n \ -1/2 



l-xf{v) /. \ -1/2 

= 2V2 l^-v^ - V,{x)j dx . 

By the change of coordinates z = x/xf\v), we obtain that 



Hence, we get 

r,(o) : 

Therefore, we have that 



r,(0) = 2vr, ^(0) = 0, 1^(0) = -^, 0(0) = 



Te{v) = 2TT-^iv^±0{v^) . 



□ 



Proposition 7.3. Let V{x) be a a perturbed quartic potential. The period 
map T : {—vo,vq) — > R satisfies the following estimate 

T{v) = 2TT-^iv^±0{\v\^) . 

Proof. By Lemma 16.31 there is f > small enough, T > and C > with 
the following properties: for all \v\ < vq, we have that 2Ti{v) < T, and for 
all < t < r, we have that 

\x{v\t) - X(>{v]t)\ < Cv^ . (37) 

By Lemma l5.8|. there is -ftT > such that, for all |v| < uq, we have that 

xe{v;2Cv^) > 2Cv'^ - K\vf (38) 

Xi{v;-2Cv^) < -2Cv'^ + K\vf. (39) 

Combining (1370 and ()38p . we obtain that 

x{v;Ti{v) + 2Cv^) > xe{v;Te{v) + 2Cv'^) - Cv"^ 

> xe{v;+2Cv^) -Cv^ 

> Cv^ - K\vf . (40) 

Combining (f37l) and ([39l) . we get 

X {v;Ti{v) - 2Cv^) < xe{v;Te{v) - 2Cv^) + Cv'^ 

< xe{v;-2Cv^) + Cv^ 

< -Cv'^ + Klvf. (41) 

Take vi < vq such that Cv'^ > K\v\^. Combining (gOl) and dH]), for all 
\v\ < vi we obtain that 

x {v; Ti{v) - 2Cv^) < < x {v; Ti,{v) + 2Cv^) . 

Therefore, by continuity of x [v] t), we have that 

Ti{v) - 2Cv^ < T{v) < Ti{v) + 2Cv^ . 
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Hence, by Lemma 17.21 we get 



Svr 

'.v^ ± U [\v\^ ) . 

□ 



T{v) =2tt-^£v'^ ±0{\v\^) . 



8. The scaling parameter 

In this section, we will use the period map T{v) to determine estimates 
for the {n,t)-scaling velocity ^n,t with relevance for the determination of 
the (n, t)-scaling parameter Tn^t introduced in ([8]). Furthermore, we will 
compute some estimates for x{Tn,tV',t) which we be useful for the study of 
the n-renormalized trajectories Xn{v;t) = (— l)"r~^x {Tn,tV',mr — it). 

The {n,t)-scaling velocity jn.t ^ has the property that there exists 
Ki > such that 

xijn,t;nTT - it) =0 

for every < t < Kin. Furthermore, we have that x{'yn~.i^t~£Tri nn — it) = 0. 

Proposition 8.1. Let V{x) be a perturbed quartic potential and i € {—1, 1} 
be equal to the sign ofV^'^\0). There exists K > such that for every n>l 
and for every t £ [0,Kn] the {n,t)-scaling velocity ^n,t > given by 

2if 

T{jn,t) = 27T (42) 

n 

is well-defined. Furthermore, the following estimates are satisfied 

3/2\ 



(43) 



Tl 

-^T ijn,tv) = nn -it- Re^^^t ± O . , . , 

where Ri^v^t = —^t (l — v'^) . In particular, we have thatj^.t G O (^{\t\/n)^^'^ 

Proof. For simplicity of notation, we will denote 'yn,t by 7 throughout the 
proof. By Proposition 17.31 there is vq = vo{V) > such that, for all \v\ < vq, 
the period map T{v) satisfies 

T{v) =2Tr-^iv'^ ±0{\v\^) . (44) 



Hence, we have that 



r(^)-2vr = -^^72±0(|7|') . (45) 



By (IM]), we get 
Combining (j46p and (|45p . we obtain 



2it 

T{j)-27T = -— . (46) 
n 
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Therefore, there is K > such that 7 = jn,t is weh-defined, for every 
< t < Kn, and 7^ e O (^). Furthermore, 



7^ e 



Hence, by (IMj) . we get 

r(7)-r(7t;) G 



8t 

37rn 



±0 



3/2^ 



-^^(7^-7V)±0 





3/2\ 


n 





C -2£-(l-v2)±0 



n 



t 

n 
3/2 



3/2^ 



Thus, by (gS]), we have that 
nvr — It 



2^(7) 



e -T (71;) - (1 - ± O 



1/2 



Remark Note that inequahties (j43p can be restated as 



7„/ e r„/ ± o 



3/2^ 



□ 



(47) 



Lemma 8.2. Let V{x) be a perturbed quartic potential and i £ {—1,1} be 
equal to the sign of ^(^^(O). There exists K > such that for every n > 1 
and for every t € [0, Xn] the {n,t)-scaling parameter Tn^t satisfies 



n7r-£t = -T {Tn,t v) + Re^^^t ± O 



n 



1/2 



where Ri^v,t = —^t (l — v^) . 

Proof. For simphcity of notation, we will denote 'yn,t by 7 and Tn^t by F 
throughout the proof. By Lemma l8. 11 we get that 



mr — it 



n 



From Proposition 17.31 we obtain that 



^^7V±0(7^) . 



(48) 



(49) 



Putting together identities (|47|) and (j49|) and Proposition [731 we get 
T{^v) G 27r - ^£F V ± O (f3) 

c T{rv)±o{r^) . 



(50) 



26 C. A. A. DE CARVALHO, M. M. PEIXOTO, D. PINHEIRO, AND A. A. PINTO 

Combining 1^ and ([SOD and recalling that T e O {{\t\/n)'^/'^) , we obtain 

■|t|3/2\ 



mr-it = (Tv) - it (l - v'^) ± O . , ,^ , , 
2 V JT- ' / 

completing the proof. □ 

Lemma 8.3. Let V{x) be a perturbed quartic potential. There exists K > 
such that the solution x (Tn^tv, t) of the Euler-Lagrange equation associated 
to the potential V{x) satisfies the following estimate 

X {rn,tv; t) e rn,tv sin (t) ± o (|r„,t , 

for allve[-l, 1] and for all t E [0, Kn^^^] . 

Proof. For simplicity of notation, we will denote Tn^t by T trough out the 
proof. Let x > be such that t = xT (Tv), and denote by [x] the fractional 
part of x. By Proposition 18. H we have that x = 0{t) and 

t = {x - [x\)T {Tv) + [x]T {Tv) 

G 2t: {x - [x]) + [x]T {T v) ± O {tT'^) . 

Hence, we get 

[x\T {Tv) et-2Trn±0 {tT^) . (51) 
By Proposition 16.41 we have that 

x{Tv;t) = X (Tv; {x - [x])T (Tv) + [x]T {Tv)) 
= X {Tv; [x]T (Tv)) 
£ Tv sin {[x]T{Tv))±0{T^) . 
Hence, by (|51|) . we obtain that 

x{Tv;t) = Tvsm{[x]T{Tv))±0{T^) 

e Tvsin {t - 27r(x - [x]) ± O (tP^)) ± O (P^) 
C Tv sin {t±0 (tp2) ) ± O (P^) 
C Tvsm{t) ±0 {tf"^) . 

Choosing K > small enough, for all t G [0, Kn^^^] , we have that O {tT^) < 
O (t5/2/n3/2) < 1. □ 

9. Asymptotic universality of trajectories 

In this section we combine the results of the previous sections to prove 
Proposition 14.21 the n-renormalized trajectories Xn{v;t) converge to 
universal asymptotic trajectories X£{v;t) introduced in section [3l 



Proof of Proposition For simplicity of notation, we will denote Tn.t by 
P trough out the proof. By Lemma 18.21 we have that 

X {Tv; HIT -lt)=x {Tv; {Tv) + Ri^t,v ± O {q^n 

Since T {Tv) is the period of x (Tv; t), we get 

X {Tv; nvr - £t) = x (Tv; ^ T (Tv) + Ri^^^t ± O (P^n 
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where [n/2] denotes the fractional part of n/2. By Proposition 17.31 we have 
that T {Tv) =2tt±0 (F^) . Hence, we have that 

x{Tv;n-K - it) = x{tv;2-k ^ + ± O (F^ + T^m)) . 

By Proposition 18.11 we have that R^^v.t G 0{t). Hence, by Lemma [8^ we 
get 

x{Tv;mT-H) G (-l)"rt; sin (i?^ ^ j ± O (P^n + r2))±0 (P^ (p^ + F^n + t)) . 

(52) 

By Proposition ([8]), we have F G O (^(l^l/"-)^^^^ • Hence, we fix X > small 

enough so that O (P^n + F^) < 1 for every t G [0,/Cni/3]. Thus, using the 
trigonometric formulas, we obtain that 

sin {Ri^^^t ± O (F^n + F^)) G sin (i?£,„,t) ± O (P^ + F^n) . (53) 

Combining (|52|) and (|53|) . we obtain 

X {Tv; rnr - It) G (-l)"Ff sin (i?^,^,t) ± O (P^ + P^t + F^n) . 

Therefore, we get 

{-lYT-^x{Tv]n^ - tt) G ?;sin(i?^,^,t) ±0 (F^n + F^ + p2t) 

C Xi{v;t)±0{T'^ + TH + T^n) . 

Therefore, we obtain that 

( \t\ ItP 

(-1)"P-^X (F?;; nvr - H) G X^(t;; t) ± O ^ + ^ + ' ' ' 



n n n 



1/2 



□ 



The next result follows easily from Proposition 14.21 It is worth to point 
out that it is equivalent to Theorem 13.11 



Corollary 9.1. Let V{x) he a perturbed quartic potential and let K > be 
as in Proposition \4.^ For every < e < 1/3, we have that 

||x„(z;;t) -X^(t;;i)||po([_i,i]x[o,i^ni/3-.],R) < O (n'^'^^^ , 

where £ is equal to the sign ofV^'^\0). 

10. Conclusions 

We have studied the dynamics of a family of mechanical systems that in- 
cludes the pendulum at small neighbourhoods of an elliptic equilibrium and 
characterized such dynamical behaviour through a renormalization scheme. 
We have proved that the asymptotic limit of the renormalization scheme in- 
troduced in this paper is universal: it is the same for all the elements in the 
considered class of mechanical systems. As a consequence we have obtained 
an universal asymptotic focal decomposition for this family of mechanical 
systems. 

We believe that the existence of an universal asymptotic focal decompo- 
sition might be useful not only on the theory of boundary value problems of 
ordinary differential equations but also on several distinct fields of the phys- 
ical sciences such as quantum statistical mechanics, optics, general relativity 
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and even tsunami formation. Our belief on such applications is based on the 
relevance that the concept of focal decomposition may have on the study of 
caustic formation by focusing wavefronts, of such significance to those fields. 
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